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Combining first-principles calculations with a technique for many-body problems, we in- 
vestigate properties of the transition metal oxide Sr 2 V04 from the microscopic point of 
view. By using the local density approximation (LDA), the high-energy band structure is 
obtained, while screened Coulomb interactions are derived from the constrained LDA and 
the GW method. The renormalization of the kinetic energy is determined from the GW 
method. By these downfolding procedures, an effective Hamiltonian at low energies is de- 
rived. Applying the path integral renormalization group method to this Hamiltonian, we 
obtain ground state properties such as the magnetic and orbital orders. Obtained results 
are consistent with experiments within available data. We find that Sr2V04 is close to the 
metal-insulator transition. Furthermore, because of the coexistence and competition of fer- 
romagnetic and antiferromgnetic exchange interactions in this system, an antiferromagnetic 
and orbital-ordered state with a nontrivial and large unit cell structure is predicted in the 
ground state. The calculated optical conductivity shows characteristic shoulder structure in 
agreement with the experimental results. This suggests an orbital selective reduction of the 
Mott gap. 

KEYWORDS: Sr2V04, local density approximation, first-principles calculation, path-integral 
renormalization group, spin-orbital coupled order, orbital degeneracy 



1. Introduction 

The discovery of high temperature superconductivity in cuprates has stimulated interest 
to strongly correlated electron systems, such as transition metal oxides and heavy fermion 
compounds. These materials exhibit a variety of magnetic, dielectric, and transport properties, 
which are caused by the strong Coulomb interactions between conduction electrons. 1 

Sr2V04 is one of the typical transition metal oxides isostructural to the parent compound 
of the high-T c superconducting material, La2Cu04, which has K^NiF^type structure (shown 
in Fig. 1(a)). Namely, this compound has a layered perovskite structure with a strong two- 
dimensional anisotropy. The similarity (or duality) of the electronic configuration between V 4+ 
(one 3d electron) in Sr2V04 and Cu 2+ (one 3d hole) in La2Cu04 should be noted. However, 
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this duality is not complete. The crystal-field effect splits 3d orbitals of V 4+ into e g and t2 g 
orbitals (see Fig. 1(b)). Furthermore the strong two-dimensionality splits t2 g orbitals into the 
xy orbital and degenerate yz and zx orbitals where the orbital degeneracy of d 1 electron 
remains between yz and zx orbitals, in contrast to the orbitally non-degenerate structure of 
the cuprates, where d x i_ y i orbital is isolated near the Fermi level. In the cuprates, in addition, 
the oxygen p orbital is strongly hybridized while this hybridization is small in Sr2V04. The 
crystal field splitting between xy orbital and others (yz and zx) are also rather small (~ 0.08eV 
in the LDA calculation) in Sr2V04 and the xy orbital is also expected to play some role in the 
ground state. This comparison suggests that studies on Sr2V04 may contribute in revealing 
the role of orbital degeneracy in comparison to the cuprates. It may also provide us with 
insights on the physics and mechanism of the high-T c superconductors. 

Motivated by the analogy with La2Cu04, there exist several experimental studies on 
Sr2VC>4. 2 ~ 7 Recently thin film on LaSrA104 became available and investigated. ' 8 The elec- 
tronic resistivity shows semiconducting behavior, whereas the resistivity does not diverge at 
low temperature. 5 ' 6 The activation energy observed in the temperature dependence of the 
resistivity lies in the range of 50-90 meV 5 and 70 meV, 8 which indicates that the Sr2V04 
is a Mott insulator with a very small Mott gap. In fact metallic phase easily appears by 
La doping. ' 8 The magnetic susceptibility follows the Curie- Weiss law at temperature above 
100K. With decreasing temperature, cusps appear at 45K 5 and at 10K. 2 Although the value 
of Neel temperature strongly depends on samples, the ground state appears to possess an 
antiferromagnetic (AF) order. Because these experimental results were obtained mostly from 
powder samples except for thin film results, ' 8 it is hard to draw conclusions on the intrinsic 
properties. As far as we know, studies on orbital orders are not available in the literature. 
The optical conductivity measured for thin films show a characteristic shoulder structure, 7 
which requires a careful consideration because it is not easily explained from the Hartree-Fock 
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Fig. 1. (a) Crystal structure of Sr 2 V04 and (b) 3d levels. 



2/38 



J. Phys. Soc. Jpn. Full Paper 

picture. 

To get insight into physics of Sr2V04 and predict intrinsic properties of this compound 
in an accurate computational framework, we apply a newly developed computational method 
proposed as a hybrid first-principles approach combined with the low-energy solver for the 
low-energy effective model. 9 For the low-energy solver, we employ the path-integral renor- 
malization group (PIRG) method. 10-14 As we will remark later, more conventional theoret- 
ical approaches predict rather controversial results: The local density approximation (LDA) 
concludes paramagnetic metal while the Hartree-Fock approximation (HFA) suggests the fer- 
romagnetic insulator. In the present paper, we show that these two results are not accurate 
and the careful considerations on the strong electron correlation effect and quantum fluc- 
tuations are crucially important in analyzing the properties of this compound. The present 
first-principles PIRG calculation predicts that this compound is very close to the Mott tran- 
sition in the insulating side with a tiny Mott gap. In addition, it predicts that the spin and 
orbital are strongly coupled with frustration effects, leading to a nontrivial and long-period 
ordering structure of spins and orbitals. We also show that the optical conductivity is consis- 
tent with the experimental results and is interpreted by the orbital selective reduction of the 
Mott gap. A part of the present results has already been reported in a letter. 9 

This paper is organized as follows. In the next section, we derive the Hubbard type effective 
Hamiltonian from the results of the first-principles calculation. In §3 and §4, we refer to the 
low-energy effective Hamiltonian and the conventional approaches. In §5, the formalism of 
PIRG method and its application are briefly described. We show results by HFA and PIRG 
and discuss physical properties of Sr2V04 in §6. A brief summary and discussions are given 
in §7. 

2. High-Energy Electronic Structure and Downfolding 

Since electrons far from the Fermi level are not usually affected by strong fluctuations 
arising from correlation effects, the LDA and GW approximation offers reasonable and suitable 
methods in this regime. We perform elimination of the high-energy degrees of freedom by 
a downfolding procedure by the density functional approach. We then derive the effective 
low energy Hamiltonian, 15-17 which can be treated by the model Hamiltonian approach. We 
employ a downfolding procedure, in which the LDA based on the LMTO basis functions is 
combined with the GW approach and the constrained LDA method to derive the screened 
Coulomb interaction and the renormalized band dispersion. 

Since global band structure obtained by LDA calculations shows that 3d t2 g orbitals are 
close to the Fermi level and are rather isolated from other orbitals (Fig. 2), it allows to 
construct the effective low-energy Hamiltonian by extracting only the degrees of freedom for 
the V 3d t2 g in Sr2V04. Since the crystal field splitting between xy orbital and other orbitals is 
quite small (~ 0.08eV in LDA calculation), xy orbital should not be neglected. We formulate 
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Fig. 2. (color online) LDA density of states. The energy is measured form the Fermi level (0 ev). 



the downfolding procedure in detail in this section. 

2.1 Kinetic- Energy Part 

The kinetic-energy part of the model Hamiltonian has been derived using the downfolding 
method, 15,17 starting from the electronic structure in LDA. 

For practical purposes we use the linear-mumn-tin-orbital (LMTO) method supplemented 
with the atomic spheres approximation (ASA). 18 The LMTO basis functions, {|x)} (the- so- 
called muffin-tin orbitals - MTOs) have many similarities with the orthogonalized atomic 
orbitals. They are constructed from solutions of Kohn-Sham (KS) equations inside atomic 
spheres, (j), calculated at certain energies E v (typically, the center of gravity of occupied band 
or of the entire band), and their energy derivatives <fi. At the atomic spheres boundaries, <p and 
cf) match continuously and differentiably onto certain envelope function (typically, irregular 
solutions of Laplace equation, which rapidly decay in the real space). Since (cp\(j))=0, the basis 
is orthogonal and {|x)} can be regarded as the Wannier functions for the full KS Hamiltonian 
Hks- 

In order to describe properly the electronic structure of Sr2V04 in the valent part of the 
spectrum, we used 67 MTOs (including those associated with the empty spheres, which have 
been added in order to improve the atomic spheres approximation for the loosely packed lay- 
ered perovskite structure - see Table I). The KS Hamiltonian in the basis of these MTOs will 
be denoted by H={x\Hks\x) ■ Corresponding densities of states and the energy dispersion of 
the t2 g bands are shown in Figs. 2 and 3, respectively. We note an excellent agreement with cal- 
culations employing more accurate full-potential linearized-augmented-plane-wave (FLAPW) 
method. 19 

What we want to do next is to describe the behavior of three t2 g bands located near the 
Fermi level by certain tight-binding (TB) Hamiltonian h, which, contrary to H , is formulated 
in basis of only three (unknown yet) Wannier orbitals centered at each V site. 
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Table I. Atomic positions (units of a=3.837Aand c=12.576A), atomic radii (in A) and basis functions 
included in LMTO calculations of Sr2V04. The symbol 'Em' stands for the empty spheres. The 
crystal structure is shown in Fig. 1. 



type of atom position atomic radius LMTO basis number of atoms 



Sr 


(0,0,0.354) 


1.924 




2 


V 


(0,0,0) 


1.383 


AsApM 


1 


Oi 


(0.5,0,0) 


0.988 


2s2p 


2 


o 2 


(0,0,0.158) 


1.012 


2s2p 


2 


Emi 


(0.5,0,-0.25) 


0.648 


Is 


2 


E1112 


(0.313,0,-0.208) 


0.437 


Is 


8 



By denoting the subspace spanned by three MTOs of t2 g type as {\d)} and the rest of 
MTOs as {|r)}, such that {|x)} = {|d)}©{|r)}, the equations for eigenvalues and eigenfunctions 
of the KS Hamiltonian H can be rearranged identically as 

{H dd -uj)\d) +H dr \r) = 0, (1) 

H rd \d) + {H rr - u)\r) = 0. (2) 

By eliminating \r) from eq. (2) and substituting it in eq. (1) one obtains the effective in- 
dependent Hamiltonian in the basis of d-states 

H d i(u) = H dd - H dr (H rr - uj)- l H rd 

and the "overlap" matrix 

S{u) = 1 + H dr (H rr - cuy 2 H rd , 
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Fig. 3. (color online) The band structure of V 3d t^g orbitals for Si^VOzt. Solid lines represent the 
result of LMTO, while dashed-dotted lines stand for our tight binding approximation. 
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such that (g?|iS'|cZ)=1. Then, the required TB Hamiltonian, h is obtained after the orthonor- 
malization of the vectors \d)^\d)=S l / 2 \d) and fixing the energy ui in the center of gravity of 
t 2g band (uq): 

h = S-V^Hg^S-Wiuo). (3) 

The downfolding method is nearly perfect and well reproduces the behavior of three t2 g 
bands in the reciprocal (k) space (Fig. 3). The obtained Hamiltonian is then Fourier trans- 
formed into the real space. The site-diagonal elements of A=||/i^ m '||, obtained after such a 
transformation, describe the crystal-field (CF) splitting, whereas the off-diagonal elements 
stand for transfer integrals. Here i and j stand for atomic sites while m,m' represent the 
band indices. 



2.2 Wannier Functions 

The Wannier functions for the isolated t 2g band can be reconstructed from the parameters 
of the kinetic energy obtained in the downfolding method. 17 In order to do so, we solve an 
inverse problem and search for the functions {H 7 /™}, which after applying to the KS Hamil- 
tonian in the real space, generate the matrix elements h™ m ' ={W™\Hks\ W™'). The method 
has been described in details in the previous publication. 17 The basic idea is to search H 7 ™ in 
the form of MTO: 18 

j,m' 

with subsequent orthonormalization and orthogonalization to other bands. 

The Wannier functions obtained in such a way are shown in Fig. 4, and their extension 
in the real space is illustrated in Fig. 5. (color online) Since the t 2g band is an antibonding 
combination of the atomic V(3d-t2 g ) and 0{2p) orbitals, the Wannier functions have nodes 
located between vanadium and oxygen sites. The xy-orbitals appear to be more localized. 
This is not surprising because in order to describe the distribution of the electron density 
in the xy-plane, it is sufficient to have compact Wannier orbitals, which are confined mainly 
within the central V site and the neighboring oxygen sites. Then, the full space within the 
plane can be spanned by such Wannier orbitals, periodically repeated in the xy-plane. On 
the other hand, the yz and zx orbitals should account for the penetration of the electron 
density in rather extended interplane region, including the areas of weak and nearly constant 
densities associated with the empty spaces of the layered perovskite structure (Emi in Fig. 4). 
Therefore, the yz and zx orbitals should have a long tails spreading out of the xy-plane. 
Another measure of the localization of the Wannier functions is the expectation values of the 
square of the position operator (r 2 ), 20 which yields 1.53 and 3.18 A 2 for the xy and yz (zx) 
orbitals, respectively. 
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Fig. 4. (color online) Contour plot of Wannier functions of the xy (left) and yz (right) types for 
Sr2V04. The solid and dashed lines correspond to the positive and negative values of the Wannier 
functions. Around each site, the Wannier function increases/decreases with the step 0.04 from the 
values indicated on the graph. Atomic positions are shown by the following symbols: o (V), x (O), 
and + (Emi). 
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Fig. 5. Spacial extension of Wannier functions for S^VO^: the number of electrons accumulated 
around the central V site after adding every new sphere of neighboring atomic sites. 

2.3 Effective Coulomb Interactions 

The matrix elements of the effective Coulomb interactions in the t2 g band are defined as 
the energy cost for moving an electron between two Wannier orbitals, W™ and Wj 71 , which 
have been initially populated by tTj m and xij m i electrons: 

A™ m = E [tli m + l,% m ' - l] - E [tli m ,%- m >] • 

For v£j, A% m have a meaning of on-site Coulomb interactions screened by intersite in- 
teractions, while for i=j it stands for the intraatomic exchange and nonsphericity effects, 
responsible for the Hund's first and second rule, respectively. 

These matrix elements can be calculated in two steps, employing a number of approxima- 
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tions, which have been discussed in details in the previous paper. 17 

1. First, we perform the standard constraint-LDA (c-LDA) calculations, where all matrix ele- 
ments of hybridization involving the atomic V(3d)-states have been artificially switched off. 21 
This part takes into account the screening of on-site Coulomb interactions caused by the relax- 
ation of the 3d-atomic basis functions and the redistribution of the rest of the charge density. 
Particularly, the removal of the 3d electron creates at the same site an excessive charge of 
the 4sp-type, which compensate the lack of the 3d charge by nearly 50%. The values of the 
on-site Coulomb (u) obtained in such an approach vary from 9.7 till 10.4 eV, depending on 
the distance between sites i and j. They are considerably smaller than the bare Coulomb in- 
teraction between 3d electrons (~22 eV), meaning that the screening of Coulomb interactions 
caused by the relaxation of the basis functions as well as the redistribution of non-3<i electrons 
is very efficient and reduces u by more than 50%. The exchange interaction (j) is about 1 eV. 
Using only u and j, one can restore the full 5x5x5x5 matrix u=||u mm / m » m /» || of Coulomb 
interactions between atomic 3d-orbitals, as it is typically done in the LDA+C7 method. 23 

2. Then, we switch on the hybridization and evaluate the screening caused by the change of 
this hybridization between the atomic V(3d)-orbitals and the rest of the basis states in the 
random-phase approximation (RPA): 

U{uj) = [l - uP(w)l 1 u. (4) 

This scheme implies that different channels of screening can be included consecutively. Namely, 
the it-matrix derived from c-LDA is used as the bare Coulomb interaction in the Dyson 
equation (4), and the 5x5x5x5 polarization matrix P=\\P mm ' m "m"'\\ describes solely the 
effects of hybridization of the V(3eZ) states with the 0(2p) and Sr(4c2) states, which leads to 
the formation of the distinct V(t2 g ) band in Fig. 2. The matrix elements of P are given by 

p f,A_ 9 W ( n "k ~ tWk04m'k'^ , "k4i"nk d m'"n , k' 

nk n'k' 

where {e n k} and {n n k} are the LDA eigenvalues and occupation numbers, respectively, for 
the band n and momentum k in the first Brillouin zone, and d m „k=(Xm|V'nk) is the projection 
of the LDA eigenstate ^„k onto the atomic 3d orbital % m . 

Results of these calculations are shown in Fig. 6. In this case one can easily separate differ- 
ent contributions to the RPA screening caused by the interband 0(2p)^Y(e g ), 0(2p)^V(t2 S ), 
and V(t2 9 )^V(e 9 ) transitions in the polarization function (5). Note that the screening caused 
by the V(<2 S )— ^(£23) transitions in the polarization function is not included here because all 
Coulomb interactions within £23 band will be treated rigorously using the PIRG method, be- 
yond RPA. Since the V-states of the t2 g and e g symmetry are practically not mixed by transfer 
interactions in the two-dimensional perovskite lattice, the screening of on-site Coulomb inter- 
actions associated with the V(t2 9 )^V(e 3 ) transitions appears to be weak. On the other hand, 
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Fig. 6. Screened intra-orbital Coulomb interaction in constraint-LDA and RPA. The partial contri- 
butions show the screening caused by separate intcrband transitions in the polarization function 
(5). 

the screening caused by the 0(2p)^V(e 9 ) and 0(2p)^V(t2 g ) transitions is very efficient, so 
that the characteristic parameters of inter-orbital Coulomb interactions are strongly reduced 
in comparison with c-LDA, typically till 2.6-2.8 eV. 

This, however, implies a strong w-dependence of the screened Coulomb interaction, as it 
immediately follows from the Kramers-Kronig transformation in RPA: 24 

ReU(u)=u VI dw — 

vr J lo 2 -uj' 2 

Indeed, the change of ReU(u>) at u>=0 is directly related with the spectral weight of |ImC/(w)| 
at finite u (Fig. 7). This spectral weight can contribute to the renormalization of the low- 




1 2 3 4 5 
0>(eV) 

Fig. 7. (color online) Frequency-dependence of intra-orbital Coulomb interaction in RPA. The solid 
(dashed) lines represent real (imaginary) parts. The inset shows high-frequency part of U(oj). 
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energy part of the spectrum through the self-energy effect. The latter can be evaluated in 
the GW approach, 24-26 where the self-energy is given by the convolution of U(ui) with the 
one-particle Green function G(u) for the tig band: 

t(u) = — [ duj'G{uj + oj')U{J). 
2vr J 

Then, the renormalization factor is 

Z= [1-aReS/^Uo]- 1 , 

and the numerical estimates yield Z«0.8. The momentum dependence ot Z is negligibly small. 
The imaginary part of E in the low-frequency region is also small and can be neglected. 

A new aspect of the RPA screening in the tetragonal compound is the anisotropy of 
screened Coulomb interaction. Since the (xy) and (yz,zx) ti g levels belong to different repre- 
sentations of the point symmetry group, the corresponding matrix elements of the screened 
Coulomb interaction U can be also different. For the static interactions this anisotropy is of 
the order of 0.2 eV, where the matrix elements between more extended yz and zx orbitals are 
better screened. 

The intersite Coulomb interaction can be also derived and the effective Hamiltonian with 
the intersite Coulomb term may also be solved by our low-energy solver. However, in this 
paper, we do not consider this effect because the intersite Coulomb interaction is substantially 
smaller. For example, the value of intersite Coulomb interactions V derived from the standard 
c-LDA approach (part 1 of the above procedure) for the cubic SrV03 is about 1.2 eV, which 
is not negligible. However, by taking into account the relaxation effects associated with the 
V(t2#)^V(e s ) transition, which can be also treated in the framework of c-LDA, the intersite 
Coulomb interaction V is reduced drastically from 1.2 till 0.3 eV. 17 Furthermore, we expect 
this value to be further reduced by taking into account the remaining 0(2p)^V(t2 g ) and 
V(t2#)^V(e 9 ) transitions in the polarization functions. 

The effective low-energy Hamiltonian is now derived from the Coulomb interaction ob- 
tained in §2.3 and the transfer. The transfer is given by the downfolded tight binding parameter 
/i™' m remormalized by Z, where the self-energy effect reduces the band width. Namely, we 
take the transfer t™f = Zh™f. 

3. Low-Energy Effective Model 



The effective low-energy Hamiltonian for the tig orbital obtained by the downfolding 
described in §2 has a form of the extended Hubbard model and is given by 

H cS = H kin + H mt (6) 

tf k in = E E *r' c L^w. (?) 

<i,j> m,m' ,a 
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-> X 



Fig. 8. Transfer directions of 3d ti g electrons on V sites in Sr2V04; Open circles denote the lattice 
sites of V. Nearest neighbor sites (N.N), 2nd neighbor sites (2.N.), 3rd neighbor sites (3.N.), and 
z direction (z) transfers are depicted. 
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(8) 
(9) 



where cJ mcr (Q mcr ) is the creation (annihilation) operator of a conduction electron with spin 
a = (1M)j an< i orbital m=(xy,yz, zx: ti g orbitals) at site i. t"* m ' is the hopping matrix 



and rij 



c lma c imcr- JJ m (K m ' m ') is the intra-orbital (inter-orbital) Coulomb interaction 
and J m ' m> is the exchange interaction, which gives rise to the Hund's rule coupling and pair 
hopping. 

For Sr2V04, estimate of each parameter obtained in the procedure explained in §2 are 
given as U xy = 2.77 eV, U yz = U zx = 2.58 eV, K xy ' yz = K zx ^ = 1.35 eV, K yz ' zx = 1.28 
eV, J xy > yz = J zx ' x v = 0.65 eV, and jy z > zx = 0.64 eV. These values approximately hold the 
relation U ~ K+2J, which shows that the effective low-energy Hamiltonian satisfies rotational 
invariance in spin and orbital spaces. 

The transfer direction is shown in Fig. 8. By considering transfers until third neighbor sites, 
tight-binding approximation can reproduce the LDA band structure (Fig. 3). Off-diagonal 
hoppings for orbital indices become zero or quite small, so that we show the diagonal hopping 
amplitude in Table II. However, note that all components of hopping matrices within third 
neighbor sites are considered in our calculations. 

Because the 2D character of Sr2V04 is strong, in the present study, we consider 2D 
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Table II. Diagonal components of hopping matrices for orbital indices, where the energy is in eV. 





N.N. 


2.N. 


3.N. 


z 


j.xy,xy 


-0.218 


-0.075 


0.000 


0.000 




-0.045 


-0.010 


0.000 


-0.016 


j.ZX,ZX 


-0.194 


-0.010 


0.018 


-0.016 



N = L x L lattice first and then N = L x L x L z lattices to understand three-dimensional 
effects. 

4. Conventional Approach 

Before going into detailed results obtained from the PIRG method, we first discuss results 
obtained from more conventional approach. The density of states (DOS) of Sr2V04 calculated 
by the LDA is shown in Fig. 2. A sharp peak around the Fermi level appears in the DOS, which 
indicates that the ground state is a paramagnetic metal. A similar result has been obtained 
by Singh, et a/. 19 These results show that the LDA result contradicts the experimental results 
for Sr2V04. On the other hand, as we will see later, the Hartree-Fock approximation for the 
effective Hamiltonian, also underlying the commonly used LDA+U approach, 23 indicates that 
the ground state is a complete ferromagnetic insulator. These are again not supported by the 
experimental results. Since LDA and HFA offer completely different phases, it is desired to 
apply a more reliable method suited for strongly correlated electron systems. 

It is widely accepted that LDA is insufficient for strongly correlated electron systems be- 
cause of neglecting spatial as well as dynamical fluctuations arising from strong correlation and 
quantum effects. In order to go beyond LDA, so-called LDA+U 22,23 ' 27 and LDA+DMFT 28-31 
methods were proposed. In LDA+U, the Coulomb interactions are treated within HFA and 
the ground state is obtained within the framework of single Slater determinant, so that it is 
also difficult to take into account effects of electron correlations and quantum fluctuations. 
On the other hand, in the LDA+DMFT which combines LDA calculation with the dynam- 
ical mean field theory (DMFT), 32-34 it is possible to calculate dynamical properties such as 
spectral functions and the optical conductivity with dynamical fluctuations being taken into 
account. However, within DMFT, it is not easy to consider spatial fluctuations as well as 
possible nontrivial spatial order, because the Hubbard type lattice model is reduced to an 
effective impurity problem in DMFT with spatial fluctuations being ignored and its extension 
to a large unit cell is not easy either. 

5. Path-integral renormalization group method 

Recently, we proposed a new algorithm, the LDA combined with the PIRG method. 9 In 
order to solve the effective Hamiltonian eq. (6), we employ PIRG method. 10-14 The PIRG can 
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systematically include symmetry breakings and spatial fluctuations neglected in DMFT. PIRG 
can be applied to various systems which can not be treated by other numerical techniques such 
as the quantum Monte Carlo method because the sign problem does not appear explicitly. 
Furthermore a Slater determinant, given as the eigenvector of the HFA, is usually employed as 
the initial state of PIRG. Since quantum fluctuations are evaluated by a systematic expansion 
of the basis states in the renormalization process of PIRG, we can discuss correlation effects 
in comparison with results of the HFA in a systematic renormalization procedure. Therefore 
this new scheme enables us to investigate the ground state properties for strongly correlated 
electron systems, and is suitable method for evaluating physical properties from the first- 
principles. 

We briefly review PIRG method in this subsection. According to the path-integral formal- 
ism, the ground state wave function \ip g ) can be obtained from 

\1> g ) = lim exp[-rF]|0o), (10) 

r — >oo 

where \4>o) represents the initial state. This projection process is performed in the imagi- 
nary time direction r. In PIRG, the ground state is approximately represented by a linear 
combination of proper basis states, 

L 

h/> ff >«J>l&>- (ii) 

1=1 

Since basis states should be chosen to give the lowest energy within given dimensions of 
the Hilbert space, we have to optimize the coefficient u>i and relevant basis states \4>i), which 
constitute the process of our numerical renormalization. This procedure is symbolically written 

as 

|^) = exp(-ArF)|4 0) ). (12) 

The optimization of \4>i) and oji is performed in the following way: We employ Slater de- 
terminants as the basis functions. Candidates of the Slater-determinant basis states \4>i) are 
generated from the operation of exp (—Aril) with a small Ar to the previous basis |^ (0) ). To 
operate exp[— AtH], the decomposition exp [—AtH] ~ exp [— Ariikin] exp [—ATH- mt ] is used 
and the operation of exp [— ATH mt ] is taken by introducing a Stratonovich-Hubbard trans- 
formation. 10 This transformation generates more than one candidate for the next generation 
basis. Among the candidates, we employ the optimized basis which gives the lowest energy 
eigenvalue in the fixed subspace and discard other candidates. The eigenvalue is computed 
from the ground state of the generalized eigenvalue problem: 

e = 1 Jtj 3 , (13) 
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where 

[H] l3 = (frlHlh), (14) 

= (4>i\<t>j)- (15) 

We repeat the operation of exp[— AtH] and the truncation by selecting the optimized basis 
within a fixed dimension L. The former process expands the number of basis functions beyond 
L and the latter truncation reduces the number to L again. After sufficient number of these 
iterations, the energy and the basis functions converge within the fixed L. The fixed point of 
the renormalization process is then completed at the fixed L. 

The ground state obtained in this restricted Hilbert space is an approximate ground state 
of the full Hilbert space. Then we further gradually increase L and repeat the iteration to 
obtain the fixed point at each L. To estimate the correct ground state properties, one needs 
extrapolations of the number of states L to the dimension of the full Hilbert space. According 
to the general extrapolation procedure, the following relation 

(H) - (H) g oc AE (16) 

is satisfied, where (H) g is the energy expectation value for the true ground state and AE is 
the energy variance defined by, 

{Hf 

For an arbitary physical quantity A, similar relations hold as 



{A) - {A) g oc AE. (18) 

In this paper, we examine dependences of these physical quantities on the energy variances in 
detail. We increase L until the dependence converges reliably to the linear behavior and the 
linear fitting is taken in this large L region. 

We take this extrapolation as far as the linearities as eq. (16) and (18) are satisfied in 
the large L region. Following the above procedures, we can estimate the energy eigenvalue 
and other physical quantities in the ground state with a systematic extrapolation to the full 
Hilbert space. 

Recently, the quantum-number projection method was implemented in PIRG, 14 which 
restores the symmetry of the Hamiltonian in the renormalization process and improves the 
accuracy. For example, the extrapolated ground state energy is —66.879, which is within the 
statistical error of the result of quantum Monte Carlo calculation [E g = —66.866 ± 0.050) 
for the single orbital Hubbard model at half filling on a 6 x 6 square lattice with the nearest 
neighbor hopping t = 1 at the onsite Coulomb interaction [7 = 4. 

In this process, by its construction, there is no negative sign problem in any models, which 
appears in the quantum Monte Carlo method. PIRG enables us to accurately study physical 
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properties in various lattices and at arbitrary fillings. Therefore this method is suitable for the 
present purpose of investigating properties of an extended Hubbard-type models in comparison 
with existing numerical algorithms. 

Here note that in order to apply PIRG to the multi-orbital Hubbard model, we have to 
treat the projection processes of intra-, inter-orbital interactions, and the Hund's rule coupling. 
However, in order to reduce the computational time, we skip projection processes of applying 
the inter-orbital and the Hund's rule coupling terms in the Hamiltonian of exp[— tH\. This 
saving of the computational time does not seem to damage the efficiency of the projection 
process. Of course, the lowering of the energy evaluated in the generalized eigenvalue problem 
is evaluated by the correct form of the Hamiltonian including the exchange terms. 

Therefore, the lowest energy state is obtained if the renormalization projection procedure 
eq. (12) is taken until the convergence. In this text, we also note that the projection eq. (12) 
generates a local minimum of the Hamiltonian H within the fixed number of basis. It may 
be necessary to seek for the global minimum by taking as many as possible candidates of 
the initial states \<Pq). The real ground state may be obtained from the comparison of these 
local minima. This is particularly important if several phases with different symmetry are 
competing as in the present case. 

6. Results 

In this section, we show physical quantities for the effective low-energy Hamiltonian of 
Sr2VC>4 and discuss ground state properties. We first apply the usual Hartree-Fock approxi- 
mation to the effective Hamiltonian eq.(6). In the next subsection, we show results obtained 
by the PIRG and discuss the relevance to experimental results. 

Here, in order to study correlation effects systematically, we introduce the interaction 
parameter A, which scales the amplitudes of all the matrix elements (U, U', J) with the ratio of 
those interactions being fixed. The realistic value in the downfolding Hamiltonian corresponds 
to A = 1.0. Namely, we study the low-energy Hamiltonian by extending eq. (6) as 

H cS = #kin + \H int (19) 

to examine the correlation effects clearly. 

Since Sr2V04 has a layered perovskite structure with strong two-dimensional anisotropy, 
we first consider two-dimensional square lattices. Three-dimensional effects are studied in the 
later part of this paper. 

6.1 Hartree-Fock results 

The Slater determinants are usually employed in the HFA and for the basis states in PIRG 
method as well. Crucial difference of PIRG from HFA is that PIRG allows representing an 
approximate ground state by increasing number of Slater determinants while HFA represents 
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the ground state only by a single Slater determinant. In other words, PIRG is able to include 
quantum fluctuations by improving HFA results, so that we first study this Hamiltonian within 
HFA as a starting point. 

Here we show results obtained by the HFA for the two-dimensional square lattices. Since 
the amplitude of orbital off-diagonal hoppings is much smaller than that of diagonal hoppings, 
each base (\xy), \yz), and \zx)) is nearly orthogonal. In Fig. 9, we show several possible spin 
and orbital patterns in the S=0 and complete ferromagnetic (S'=max) sectors in the x-y 
plane, respectively. For the Fig. 9 (a) and (b) structures, typical order parameters are given 

by; o° = \ Yl,o\( n0 xyv) ~ ( n yza)\ an d rn° a = \\{ n ° a -\) ~ ( n «|)l a ^ odd-numbered sites to the 
x direction, where represents the expectation value of electron number with spin a 

and orbital a=xy,yz. Similarly, m e zx = \\{n e zx *) — {n e zx M, at even-numbered sites to the x 
direction. For the Fig. 9 (c) structures, o° = §|(n° j) — (n° x ^)| and o e = \\{riy Z ^) — (n e zx ^)\ 
at odd and even sites to the y direction are the relevant order parameters. For the Fig. 9 (d) 
structures, we take o° = \\{n xy ^) — (n° z ^}\ at odd sites to the x direction. 

These are obtained at least as metastable states of the HF solution. Because many of 
dominant orbitals are different from those of the nearest neighbor sites and rather orthogonal 
in all the panels of Fig. 9, the nearest-neighbor exchange coupling has primarily ferromagnetic 
sign because of the Goodenough-Kanamori rule. 

However, this ferromagnetic exchange is not absolutely dominant because the orbitals are 
nonnegligibly mixed and the nearest-neighbor orbitals are not strictly orthogonal. Further- 
more, substantial third neighbor hoppings between the same orbitals exist and causes the 
antiferromagnetic coupling between the third-neighbor sites. Due to the antiferromagnetic 
interaction between the third neighbor sites, frustration is induced and a complex ordering of 
Fig. 9 (a) appears. Fig. 9 (b) is also a possible order where the spin stripe structure appears 
to the x direction. 

Figure 10 shows the total energies per site in the interval 0.8 < A < 1.2. Although the 
energy for the order (b) is lower than that of (a) within HFA in S=0 sector at A ~ 1.0, the 
difference of each energy is very small (~ O.OleV). On the other hand (c) is the lowest-energy 
state for 5=max and becomes the absolute ground state within the Hartree-Fock approxi- 
mation. However, the energy of (d) state is close to that of (c), where the difference of the 
energy per site between (c) and (d) is about 0.005eV at A ~ 1.0. Because of the ferromag- 
netic and antiferromagnetic interactions coexisting in this system, the ordering structure of 
the spin and orbital may become complicated at the realistic parameter value of the material 
(A ~ 1.0). Therefore we can not exclude the possibility of other more complicated ordered 
state. Although it is difficult to determine the Hartree-Fock ground state, the unit cell of 
Sr2V04 becomes large and nontrivial in any case. 

By including quantum fluctuations by PIRG, the order (a) becomes ground state and the 
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order (c) becomes most stable in S=max sector in our calculations, which will be discussed 
in the next subsection. Therefore, hereafter we mainly consider the order (a) and (c) in the 
insulating phases in this subsection. 

In Fig. 11, we show the double occupancy, charge gap, and total energy per site in total 
S=0 sector (Fig.9(a)) and the complete ferromagnetic (S'=max) state (Fig.9(c)). The double 




(c) S=max 




Fig. 9. Possible spin and orbital patterns obtained as metastable states of HFA; the S—0 state ((a) 
and (b)) and the complete ferromagnetic state ((c) and (d)) when A ~ 1.0. Up (down) arrow repre- 
sents the up (down) spin. The dominant orbital and spin patterns are illustrated only schematically 
and the actual solution contains some linear combination of other orbitals. 
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Fig. 10. Total energy per site as a function of the interaction parameter A obtained by HFA for each 
spin and orbital pattern in Fig. 9 for the two-dimensional square lattices. 



occupancy D and the charge gap A c are defined as 

D = Jj^2 ^2 ( n irnan im , a ,), (20) 

i (m,a)<(m' ,a') 

A c = et°-e° k s . (21) 

The charge gap is defined as the difference between the energies of the lowest unoccupied state 
e%. s and the highest occupied state e£ s , where represents the dispersion of electrons 

When A = 0, D ~ 0.42 at 5=0 while D ~ 0.33 is obtained at 5=max. With increasing 
interaction parameter A, the jumps of double occupancies appear around A ~ 0.45 (S=0) 
and A ~ 0.80 (5=max) and accordingly each charge gap begins to open continuously. These 
results indicate that the 1st order metal-insulator transition occurs at A ~ 0.45 (5=0) and 
A ~ 0.80 (5=max) within the HFA. 

Within HFA, the complete ferromagnetic and insulating state becomes always more stable 
than the 5=0 state when A > 0.25. Therefore when A = 1.0, which is the realistic value 
estimated by the LDA, the HFA predicts the complete ferromagnetic and insulating ground 
state. This result does not seem to be consistent with the experimental results. However, the 
difference between the total energy of 5=0 and that of 5=max is the order of 10~ 2 eV because 
of the competing ferromagnetic and antiferromagnetic interactions. This system appears to 
be very sensitive to the control of electron correlations. To obtain ground state properties, we 
have to consider effects of quantum fluctuations, which are not included within HFA. 

Next, let us consider the magnetic and orbital orders. Equal-time structure factors are 
defined as 

s w = ^E^-^) ex p(^-^-^))' ( 22 ) 
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Fig. 11. (a) Double occupancy D, (b) charge gap A c , and (c) total energy per site as a function of 
the interaction parameter A obtained by HFA for the two-dimensional square lattices. The solid 
line represents S=0 state in Fig. 9 (a) and the dashed line is 5=max state in Fig. 9 (c). Note that 
the configuration Fig. 9 (a) gives the lowest energy among 5=0 states. 



T (4) = ^E< T *- T i> ex P( i 9-( r <- r i))> ( 23 ) 
where Si and T, are the spin and orbital operators as follows, 

S i = g Yj C lma a «y C irnv>, (24) 
m,cr,a' 

T% = — j= ^2 c \mij^"i,m'Ci m 'cr, (25) 



N 



where a and A are the Pauli and the Gell-Mann matrices, respectively. 

Figure 12 shows peak values of the spin and orbital structure factors by HFA when we 
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0.2 




Fig. 12. Largest spin-spin and orbital-orbital correlation functions (equal-time structure factors) by 
HFA as a function of 1/VW for the various wave vectors for the two-dimensional square lattices; 
(a) Fig. 9 (a) in total 5=0 sector at A = 0.40 (metal) and A = 1.00 (insulator), (b) Fig. 9 (c) 
in total 5=max sector at A = 0.72 (metal) and A = 1.00 (insulator). In each panel, circles and 
triangles denote the spin-spin and orbital-orbital structure factors, respectively. 

assume 5=0 and 5=max states, where wavevectors q dependences of each structure factor in 
both metallic phases are quite weak. With increasing system size, each correlation function 
obtained by the HFA monotonically decreases. All the scaled structure factors S(q)/N or 
T(q)/N of metallic phases in each spin sector become zero when N — > oo, which indicates 
metallic phases become paramagnetic (5=0) and paraorbital states (5=0 and 5=max). How- 
ever, the magnetic and orbital orders of the insulating states appear to survive when N — ► oo 
and the spin and orbital long-range order may exist in Sr2V04 within HFA of states illustrated 
in Fig. 9 (a) and (c). 

In order to investigate effects of the spin-orbit coupling, we add the spin-orbit interaction 
term H$o = CYli^i " s i to the Hamiltonian eq.(6), where ( represents the strength of the 
spin-orbit coupling. Note that we first consider z component of the spin-orbit coupling in 
this paper, for simplicity and an order parameter is defined as = —i(c^ Ba cccr — c^ Ca .CBa), 
where B and C represent yz and zx orbitals, respectively. Furthermore, by rotating each 
axis (x — ► y, y — > z, and z ^ x), we investigate the effect of x component of the spin-orbit 
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(a) 



-0.8 




Fig. 13. (a) Ground state energies per site as a function of the spin-orbit interaction parameter £ 
obtained by the HFA at A = 1.00. The open (closed) circle, inverse triangle, and square represent 
structures with the spin periodicity lxl, 2x2, and 4x2 illustrated in (b) with z (x) component 
of spin-orbit coupling. The open and closed diamond (triangle) denote (a) and (b) ((c) and (d)) 
of Fig. 9 with z component of spin-orbit coupling, respectively, (b) Structures of total angular mo- 
menta, which correspond to up and down arrows, at each wavevector q when spin-orbit interactions 
are taken into account. 

coupling. Within HFA, the ground state energy is shown as functions of £ in Fig. 13. Note that 
decoupling processes of the Hamiltonian for various ordered states are different from that of 
the Fig. 9 (a), so that each energies at ( = 0.0 do not correspond to that of Fig. 9 (a). 

Although all the energies decrease with increasing the spin-orbit coupling, those are insen- 
sitive at £ < 0.1 eV. The realistic values of spin-orbit interactions of typical transition metals 
are usually 0.01~ 0.02 eV. Even if the spin-orbit coupling is increased, the intersection of the 
energies of (a)-(d) in Fig.9 does not appear at £ < 0.2 eV, which is seen in Fig.14. Furthermore 
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Fig. 14. (color online) Ground state energies per site as a function of the spin-orbit interaction 
parameter £ obtained by the HFA at A = 1.00. The notations from (a) to (d) have one-to-one 
correspondence with those in Fig. 9. 

the energy difference per site between the ordered state in Fig. 9(a) and other ordered states 
with order parameter A so is larger than 0.15eV, so that we conclude the spin-orbit coupling 
is irrelevant to the ground state property of Sr2V04. The PIRG result shown later does not 
alter this conclusion. Hereafter we neglect effects of the spin-orbital coupling. 

The HFA results for the downfolded effective Hamiltonian of Sr2V04 predict that the 
ground state is complete ferromagnetic and insulating with the charge gap ~ 0.2 eV. These 
are not consistent with the experimental results, where the ground state appears to be an 
antiferromagnetic insulator. 

6.2 PIRG Results 

6.2.1 Metal- Insulator Boundary and Orbital-Spin Structure 

In this subsection, we show PIRG results and discuss correlation effects with quantum 
fluctuations beyond the HFA. Because of the layered structure, the coupling along the z 
direction is very small. Therefore we first consider two-dimensional lattices. Note that the 
interaction parameter A is changed at 0.08 intervals. Namely, A=0.88, 0.96, 1.04, 1.12, etc. in 
PIRG calculations, and we carry out interpolations by using obtained results at various A. 

The HFA results show that the ground state belongs to the total S'=max sector (Fig. 9 (c)), 
whereas the ordered state in Fig. 9 (b) becomes stable within the total S=0 sector. However, 
energies of other ordered states are also close to those ordered states. To estimate quantum 
fluctuation effects on this severe competition, we perform PIRG calculation starting from 
initial eigenvectors of (a)-(d) in Fig. 9 and metallic states obtained by the HFA. Even after the 
PIRG procedure, all of these states remain metastable or stable. The metallic states remain 
paramagnetic and paraorbital and ordered metallic states do not appear even after PIRG 
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Fig. 1 5 . Extrapolations of the total energy per site to zero energy variance for two-dimensional 4x4x1 
lattice at A=1.04. L is taken up to L = 192. In the upper panel, open circle, closed circle, and 
triangle represent the metal, insulator in Fig. 9 (a), and insulator in Fig. 9 (b), respectively at 
S=0 sector. On the other hand, in the lower panel, open circle, closed circle, and triangle represent 
the metal, insulator in Fig. 9 (c), and insulator in Fig. 9 (d), respectively at S'=max sector. 



calculations. In the large L region, linear fittings against the energy variance in the S=0 and 
,S=max sectors for various states are shown in Fig. 15 at A = 1.04, for example. 

Compared with results of the S=max sector, the linear fitting in the S=0 sector becomes 
satisfied only at large L since required number of basis to obtain accurate eigenstates of the 
Hamiltonian generally becomes large, which means that HFA is insufficient to describe the 
true ground state in 5=0 sector. Furthermore, this behavior is especially remarkable in the 
metallic state. It is because metallic phase tends to become unstable at large A if L is increased 
and collapses to insulating ordered states. Even at A ~ 1.0, though the metallic states remains 
metastable, the linear fitting requires larger L than that of the insulating phase. 

After this linear fittings, extrapolated values are obtained. Table III shows coefficients 
of the linear extrapolations for various states at A = 1.04. The energy holds the relation 
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Table III. Coefficients of the linear extrapolation for various states at A = 1.04 (E = sAE + E g ). 



state(5=0) 


s 


E 9 


error bar 


Metal 


13.04 


-1.105 


± 0.0025 


(a) of Fig. 9 


19.95 


-1.117 


± 0.0013 


(b) of Fig. 9 


14.58 


-1.079 


± 0.0008 


state(5=max) 


s 


E 9 


error bar 


Metal 


10.28 


-1.086 


± 0.0001 


(c) of Fig. 9 


18.64 


-1.074 


± 0.0003 


(d) of Fig. 9 


25.58 


-1.064 


± 0.0002 



E = sAE + Eg where s represents the slope of the linear extrapolation and E g is the ground 
state energy. Here note that the number of Hilbert subspace L is mainly taken up to L = 192. 

Let us discuss the total energy of metallic and insulating states for each total spin sector. 
First we consider the 5=0 insulating states where HFA result predicts that the lowest state is 
Fig. 9 (b). However, by including quantum fluctuations by PIRG, the insulating state Fig. 9 (a) 
becomes stable in comparison with Fig. 9 (b) states at A ~ 1.0. Since the spin configuration 
of Fig. 9 (b) has a stripe structure with the width 2, the ferromagnetic pattern becomes 
dominant in the x direction within HFA, so that the required number of bases to obtain the 
true ground state of Fig. 9 (b) is smaller than that of Fig. 9 (a). In fact, the energy variance 
of Fig. 9 (b) shown in Fig. 15 is small in comparison with that of Fig. 9 (a). Therefore, though 
the energies at the same L are comparable between (a) and (b), Fig. 9 (a) has lower energy in 
insulating state at 5=0 sector after the extrapolation to zero energy variance, where the total 
energy E g = —1.117 ± 0.0013 is obtained at A=1.04. Furthermore, compared with the total 
energy of the metallic state, where E g = —1.105 ± 0.0025, we can conclude that the insulating 
state Fig. 9 (a) becomes most stable in 5=0 sector. 

On the other hand, in 5=max sector, HFA result predicts that the insulating state Fig. 9 
(c) has lowest energy. In the insulating states, the total energy of Fig. 9 (c) obtained by PIRG, 
which is Eg = —1.074, becomes lower than that of Fig. 9 (d). However, compared with the 
metallic state, Fig. 9 (c) becomes a higher energy state. Therefore the metallic state becomes 
most stable in 5=max sector, where E g = —1.086 at A=1.04. 

At A ~ 1.0, in the 5=0 sector, the ground state turns out to be insulating and ordered 
as (a) in Fig. 9 after the extrapolation to zero energy variance. On the other hand, although 
Fig. 9 (c) becomes most stable in the insulating states, the ground state becomes metallic in 
the 5=max sector. Therefore, hereafter we only consider the configurations Fig. 9 (a) at total 
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Fig. 16. Total energies as a function of interaction A by PIRG calculations: total S=0 and S=max 
states for 4x4x1 lattice, where insulating phases in each spin sector correspond to Fig. 9 (a) and 
(c), respectively. 



5=0 and Fig. 9 (c) at total 5=max in the insulating states. 

Next, let us consider the metal- insulator transitions in each total spin sector. The total 
energy per site as a function of the interaction A is shown in Fig. 16 for 4x4x1 lattice. 
Although the energy of total 5=max obtained by HFA is lower than that of 5=0, the PIRG 
predicts opposite around A ~ 1.0. The energies of metallic and insulating states intersect at 
A ~ 1.0 (for the 5=0 state) and at A ~ 1.3 (for the complete ferromagnetic state), where the 
metal-insulator transitions occur. In particular, total 5=0 state is close to the metal-insulator 
transition, which indicates the realistic parameter for Sr2V04 is on the verge of the Mott 
transition. 

In order to consider the energy gain in comparison with HFA, we show the lowest total 
energies obtained by PIRG for each spin sector as a function of interaction A in Fig. 17 ((a) 
at total 5=0 and (c) at complete ferromagnetic state in Fig. 9) . The difference between the 
ground state energy of total 5=0 and lowest energy of total 5=max states is 0.03eV around 
A ~ 1.0 within PIRG results. We therefore conclude that the true ground state becomes total 
5=0 state in contrast to the result of HFA prediction. The available experimental results are 
consistent with our obtained result, where Sr2V04 shows antiferromagnetic Mott insulator 
with small gap. 

Since the candidates of configuration for spin and orbital order at 5=0 are in severe 
competition within 0.01 eV at A ~ 1.0, where the total energy of other ordered state is close 
to that of (a) in Fig. 9 (not shown), we can not exclude the possibility of other ordered states. 
However, all the candidates have large unit cells with nontrivial lattice structures. 

The obtained results are qualitatively different from the HFA results. The origin of this 
discrepancy is interpreted as follows; It is well known that the ordered state tends to be 
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Fig. 17. Lowest energies per unit cell as a function of interaction A by HFA and PIRG for total S=Q 
and complete ferromagnetic states. 

overestimated in the HFA since quantum fluctuations are not included. By taking account of 
quantum fluctuations by PIRG, the critical values of interactions for ordered states become 
large. In particular, the energy gain by quantum fluctuations is conspicuous in S=0 state than 
the ferromagnetic state. Therefore the required number of bases becomes larger in comparison 
with that of the complete ferromagnetic state to obtain an accurate estimate of the ground 
state, so that the energy gains of S=0 sector tends to become large with increasing the 
Hilbert subspace L. Therefore the single Slater determinant approximations, such as HFA 
and LDA+U approach, can not describe the antiferromagnetic states. 

On the other hand, energy gains at total 5=max by the quantum fluctuations are rather 
small because ferromagnetic configurations may be relatively well given by a single Slater 
determinant. In general, the ferromagnetic state is represented by fewer number of basis 
functions and quantum fluctuations are small. 

We conclude that Sr2V04 is close to the Mott transition and may have complicated spin 
and orbital order with Fig. 9 (a) structure for 4x4x1 lattice. Here we investigate the system 
size dependence of metallic and insulating phases illustrated in Fig. 9 (a) in S=0 sector. In 
Fig. 18, we show PIRG results of total S=0 for different lattice sizes. This result for 8x4x1 
lattice is similar to that of 4x4x 1 lattice, where insulating state for two system sizes becomes 
stable in comparison with metallic states and the metal-insulator transition occurs at A ~ 1.0. 
Therefore we believe that size dependence is small. 

Furthermore, in order to consider the three dimensional effect, we have calculated the 
total energy of some three-dimensional ordered patterns with the order (Fig. 9 (a)) on two- 
dimensional plane being fixed. However the difference of each energy is very small (within 
O.OOleV). Since the ratio of magnitudes of inter-layer hoppings and intra-layer hoppings is 
smaller than 1/10, the inter-layer correlation is quite weak. In addition, frustration is induced 
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by ferromagnetic and antiferromagnetic interactions between layers. Therefore the ground 
state has severe degeneracy, so that we can not determine the unique spin-orbital ordering 
structure in three-dimensional system within the present study. It is even conceivable that the 
interlayer correlation is frozen to short-range pattern in realistic experimental conditions. 

Here we refer to the magnetic and orbital orders. We have calculated structure factors eq. 
(22) and (23) within the PIRG for 4x4x1 and 8x4x1 lattices. Although correlation functions 
have larger error bars than the estimate of the ground state energy in PIRG and the accuracy 
of the estimate of the long-range order becomes worse, the amplitude of the long-range order 
does not have serious system size dependence and similar to the size dependence obtained in 
HFA. Therefore we believe that the magnetic and orbital orders of Fig. 9 (a) survives in the 
thermodynamic limit. 

Figure 19 shows the spin-orbital structures calculated by PIRG. (a) represents the ground 
state, which is obtained by using the HFA result (Fig. 9 (a)) as a initial state. In comparison 
with HFA structure, minor orbital component rather than the dominant orbital at each site 
becomes larger by including quantum fluctuations in PIRG, so that obtained orbitals rotate 
from each axis. On the other hand, figure 19 (b) is one of the other candidates. The energy 
difference between (a) and (b) is in the order of O.OleV, where low energy states are in severe 
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competition. 

Since the wave function can be directly calculated in PIRG, we can straightforwardly 
calculate the charge density in the realistic system after transforming to the real space repre- 
sentation of the Wannier orbitals of each " £23" electrons. Figure 20 shows the charge density 
of unit cell in Sr2V04 for each number of basis state L at S=0 state in the VO2 plane at 
z = 0, where the insulating state is indeed realized at A=1.0 (shown in Fig. 18). 

6.2.2 Optical Conductivity 

Recently the optical conductivity was measured by Matsuno et al. 7 > 8 Here we discuss the 
optical conductivity for the spin-orbital ordered state obtained by PIRG. 

Since it is difficult to evaluate dynamical quantities directly within PIRG, to calculate 
the conductivity, we complementarily use a HFA eigenstate which simulates the PIRG results 
instead of the PIRG wavefunction itself. In the PIRG results, order parameters obtained by 
HFA are quantitatively reduced because of quantum fluctuations. Therefore, though it is not 
complete, we employ the solution of HFA which reproduces the peak amplitude of the order- 
parameter structure factor of the PIRG state, by taking reduced interaction parameters. In 
order to obtain this approximate ground state, we consider the additional procedure in HFA. 
We first employ the eigenstate of HFA at A = 1.0, |3>), which is the metastable state with the 
same spin-orbital configuration as the PIRG ground state. Namely, we employ the Hartree- 
Fock state shown in Fig.9(a). Then by introducing new interaction parameters A' , A' , and 



(a) Ground State 




(b) Other candidate 




Fig. 19. Obtained spin-orbital structres when A ~ 1.0 at L=160. (a) and (b) correspond to the ground 
state and one of the other candidate, respectively. 
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Fig. 20. Contour plot of the charge density of vanadium " i 2 g" Wannier orbitals in the mirror ab 
plane consisting of V and 0, constructed for the 4x4x1 unit cell at A = 1.04 for number of basis 
at L = 160. In these contour lines, the minimal value of the charge density shown is 0.004, which 
increases with the interval 0.022. For certain oxygen sites the maximal intensity appears to be 
smaller than 0.004 and are not explicitly visible (blank spaces on the figure). 



X' zx , we newly calculate matrix elements ($\H e g\$), where the interaction part of H c g is 
replaced as, 

U X y,yz{( n xy) n yz + n xy(n yz )^j 

~l~ ^yz,zx(^{ n yz) n zx ~T" n yz( n zx)^j 

+ u 'zx,xy(( n zx)n xy + n zx (n xy f) 

~^ ^xy{u X y,yz{ n yz) + U zx ^ xy (n zx )^j n xy 

+ \z(^yz,zx( n zx) + U xy ^y Z (n xy )^jn yz 

+ ^zx(^zx,xy( n xy) + ^ yz,zxk n yz)^ n zx- (26) 

Here, U' represents inter-orbital interactions K — J or K. By diagonalizing the modified HFA 
Hamiltonian, new eigenstates \^ n ) can be obtained. By optimizing the parameters X' xy ,X' yz 
and X' zx , we seek for the Hartree-Fock ground state |^o) which reproduces relevant physical 
quantities obtained by PIRG. 

As the relevant physical quantities, we employ the charge-charge correlation functions for 
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Table IV. The largest values of charge-charge correlation functions obtained by HFA and PIRG at 
q = (n, 0) at A = 1.0 for 4x4x1 lattices. 





HFA 


PIRG 




0.45 


0.22 


Cj/z,j/z( ( l) 


0.98 


0.97 




3.1 


2.3 



each orbital m, which are defined as 

C m , m (q) = ^ Y^(rn m n jm )^- n '\ (27) 

i,j 

where rn m = n j m | + n i m ± . The largest values appear at q = (tt, 0) in all correlation functions in 
HFA and PIRG (shown in Table IV), whose results show that quantum fluctuations suppress 
each amplitude of correlation functions. 

When X' xy = 0.59, X' yz = 0.90, and X' zx = 0.62 in eq. (26), the peak values obtained by 
the modified HFA nearly coincide with the PIRG results. Therefore we regard this eigenstate 
obtained by the modified HFA as an approximate ground state. 

Then we calculate optical conductivity using this Hartree-Fock results. The optical con- 
ductivity is defined as, 

where VI/o(*m) is the ground (excited) state with the total energy Eq (Em) obtained by 
the above mentioned procedure. The current operator j a along a direction is written as 
3a = -iEt.m.m'^Mr'^L^mV ~ 4' a m> a ™*) wh ere i' x = (i x + l,i y ) and i' y = (i x ,i y + 1). 

Figure 21 shows optical conductivities at low energy regime (u < 1 eV, (a) HFA and 
(b-d) modified HFA), which consists of the contribution of excitations in the inter-orbital 
interactions K — J terms. In the higher energy region in (a), although there are other peaks 
caused by the scattering of the inter-orbital interactions (K terms, uj ~ 1.5eV) and intra- 
orbital interactions (U terms, oj ~ 3.0eV), the amplitude is relatively smaller than that of 
K — J terms. Therefore these peaks are not essential in our discussion. The K — J peak 
basically represents transitions of an electron to an orbital m on the neighboring sites where 
another electron with the parallel spin already exists at the different orbital from m. This 
process makes the energy cost of inter-orbital onsite interaction K reduced by the Hund's 
rule coupling J. The total conductivity (b) is the obtained by optimizing the parameters X' xy , 
X' yz , and X' zx , where peak positions shift to the low energy region in comparison with the HFA 
result (a). This result may approximately represent the realistic conductivity of Sr2V04. 

In order to clarify contributions from each orbital in the K — J peak of the optical conduc- 
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Fig. 21. Optical conductivity obtained by (a) HFA and (b-d) modified HFA at A = 1.0, X' xy = 0.59, 
\' yz = 0.90, and X' zx = 0.62 to simulate the PIRG result. Panels (a) and (b) are total conductivities, 
(Ttot(^) = (f xx (o;) + a yy (u)))/2 obtained by assuming the twin structure of orbital-spin order 
obtained by the configuration of Fig. 19(a) mixed with the 90 degree rotated structure. Panels 
(c) and (d) show <j xx {lo) (a yy (u)), where solid line is full conductivity and dashed, dotted, and 
dashed-dotted lines represent contributions of xy, yz, and zx orbitals. 



tivity, the transfers except for one orbital are switched off in current operator (Fig. 21 (c) and 
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(d) ). Note that the amplitude of sharp peak in a xx {to) {lo ~ 0.65 eV) is approximately twice 
as large as that in a yy (ui) because the spin-orbital structure on Fig. 19 (a) has the anisotropy 
about the direction. Since the orbital ferromagnetic structure to y direction, which consists 
of zx orbitals, is partially realized, the current along y direction is reduced. While spectra of 
the xy orbital are broad and have a diminished gap, those of the zx orbital in a xx {oj) and yz 
orbital in a yy {u) have sharp peak structures because the xy band is broader than yz and zx 
orbitals. We note that the gap of the band mainly arising from the xy orbital is significantly 
reduced as compared with that of the HFA at A = 1.0. This selective reduction of the xy gap 
makes a shoulder at uo ~ 0.40eV and a peak at to ~ 0.65eV in the optical conductivity at low 
energy region, whose overall structure is consistent with the experimental result. We note 
that it is difficult to explain the shoulder structure without the present clarification. 

There exist some disagreement of peak positions between the present result and the experi- 
mental result. In the experiment, the peak position is around 1 eV, while our result shows peak 
at smaller energy. Such discrepancy may be ascribed to the oversimplification of the PIRG 
result by the single Slater determinant in the imitated HFA. In general, even if the long range 
order disappears, the insulating gap does not always vanish. ID Hubbard model is a typical 
example. However, within the HFA, when the long range order disappears, it always become 
a metal. When we reproduce the correct order parameters by optimizing the parameters X' xy , 
X yz , and \' zx within the Hartree-Fock approximation, it ignores the gap formation through 
the genuine Mott origin without the long-range order. Then it may shift the conductivity gap 
to too low energy region, which may be the main reason of disagreement of peak positions. 

When the transfer is increased relative to the interaction from the experimental condition, 
(namely A is reduced), for example, by pressure, the present result opens an interesting possi- 
bility of orbital selective metallization, where electrons with xy component becomes metallic 
while the other component keeps the gap structure. 

Orbital selective Mott transitions (OSMT) have extensively been studied in systems with 
orbital degeneracy for d 2 system. Theoretically, dynamical mean-field theory appears to sup- 
port that OSMT takes place when the Hund exchange term is included and the orbital hy- 
bridization is absent. 35 Experimentally, (Ca,Sr)2Ru04 has been proposed to be a good candi- 
date for the occurrence of the OSMT. 36 ' 37 However, (Ca,Sr)2Ru04 has a substantial tilting of 
the octahedra (so called GdFe03-type distortions) in the perovskite structure because of the 
large ionic radius of Ru atom. This tilting yields substantial hybridization between orbitals, 
which could mask the real OSMT. In fact, theoretically it was suggested that the OSMT 
would be sensitively destroyed by the orbital hybridization. 

In Sr2V04, the tilting is absent and therefore the hybridization between orbital can be 
neglected. In fact, recent optical conductivity data 7 are successfully explained if only one 
orbital is on the verge of the Mott transition while other orbitals keep large gap, which is 
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indeed seen in our DFT-PIRG result. A crucial difference from (Ca,Sr)2RuC>4 is, however, 
that Sr2VC>4 is d 1 system, where at the metal insulator transition of xy orbital, other orbitals 
are basically empty. In other words, if the xy component of electrons is metallized, other 
component is not able to keep the commensurability condition strictly required for the Mott 
insulator. Therefore we do not expect the original sense of OSMT anyhow. However, in the 
analogy of OSMT, it would be conceivable to keep an insulating gap for the excitation of 
electrons in the yz and zx orbitals, if we expect a continuous evolution of the optical spectra. 
(Note that the d l electron in the ground state has substantial amount of yz and zx component 
in our PIRG result.) Our DFT-PIRG result, however, shows that such intermediate phase may 
be preempted by a single first-order Mott transition with concurrent vanishing of all the gaps 
of the orbitals in this realistic case. Namely, right after the transition to the metallic state, 
the PIRG result predicts that the orbital-spin order collapses and the whole gap structure 
vanishes immediately. 

Although DMFT supports the existence of OSMT for d 2 system, DMFT has a limitation 
of ignoring the spatial correlation. The present result implies that, even in the d 2 systems, 
spatial fluctuations and collapse of orbital-spin order structure may have a serious effect 
on the stability of the intermediate phase at OSMT. It would be an intriguing open issue 
whether OSMT may theoretically survive in the realistic situation, where the insulating gap is 
substantially influenced by the spatial structure of the orbital-spin order. A possible situation 
is that the collapse of one component orbital order at the metallization may destroy other 
order cooperatively, which results in the simultaneous collapse of the whole gap. 

For reference, we consider one-particle energy dispersions in our modified HFA to simulate 
PIRG result. The one-particle wave function in the momentum space is defined as, 



where n is the sublattice index (4x2 (structure) x 3 (orbitals) = 24). Ri denotes the sublattice 
site and r n is the position in the sublattice. By using the wave function in the momentum 
space, the matrix element of the modified HF Hamiltonian with k space is given by, 



Diagonalizing this 24x24 matrix for each k, one-particle energy dispersions can be determined. 
The obtained dispersions are shown in Fig. 22. The lowest 4 bands are fully occupied and other 
bands are empty, so that the small energy gap appears at ~ —0.8 eV, which is reflected in 
optical conductivity. Note that 3 bands in the fully occupied bands become quite flat because 
of the gap opening caused by the strong Coulomb interactions. 

In summarizing the PIRG results, antiferromagnetic fluctuations are larger than ferro- 



k,n) = 




(29) 




(\t,n\H\k,ri) 




(30) 
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Fig. 22. One-particle energy dispersions at A = 1.0, X' = 0.59, X' yz = 0.90, and A^ = 0.62 in the 
reduced Brillouin zone of 4 x 2 structure. Fermi level is at E F ~ — 0.8eV. 

magnetic fluctuations and these fluctuations ignored in the HFA are taken into account by 
PIRG, which leads to a larger energy lowering of the S=0 state by PIRG. The present PIRG 
result suggests that this energy lowering by the quantum fluctuations is crucial in determining 
the true ground state. We conclude the ground state becomes the insulating state with small 
gap at S=0 at A ~ 1 as seen from Fig. 16 with a nontrivial spin-orbital structure illustrated 
in Fig. 19. 

7. Summary and Discussion 

In summary, we have investigated the ground state properties of Sr2V04 by employing 
new algorithm, LDA combined with PIRG. While the high energy band structure is deter- 
mined by the LDA, the electrons close to the Fermi level are sensitive to quantum fluctuation 
effects and are studied by more accurate computational method. To bridge these two energy 
scales, we have derived the low-energy effective Hamiltonian from the high-energy LDA band 
structure by a downfolding scheme using combinations of the constrained LDA and the GW 
approximation. Since the isolated bands close to the Fermi level in Sr2V04 in LDA mainly 
consist of V 3d t2 9 orbitals, we extract these low-energy bands from the LDA bands in order 
to construct the effective Hamiltonian. Obtained effective Hamiltonian is treated by PIRG. 

Applying the HFA and PIRG schemes to the obtained effective Hamiltonian for Sr2V04, 
we have calculated the total energy and physical quantities for metallic and various insulating 
ordered states for S=0 and complete ferromagnetic states. The HFA result suggests that 
Sr2V04 becomes a complete ferromagnetic insulator with a large charge gap ~ 0.2 eV. LDA+U 
prediction has similar results to HFA. On the other hand, PIRG result shows that S=0 state 
becomes the ground state, where the large unit cell structure with complicated spin and 
orbital pattern is realized. This complicated structure is resulted from frustration effects due 
to effective long-range exchanges. Furthermore Sr2V04 is close to the Mott transition, which 
implies that the charge gap is quite small. These obtained results are consistent with available 
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experiments in contrast to the expectations of the LDA and HFA. 

While ferromagnetic states are well described by single Slater determinant approximations, 
such as LDA, HFA, et al, antiferromagnetic states with large quantum fluctuations are not. By 
considering strong quantum fluctuations by PIRG scheme beyond a single Slater determinant, 
the accurate ground state can be described. 

This material shows severe competitions. In addition to being close to the Mott transition, 
the energy differences between various complicated ordered states at S=0 sector are rather 
small. Although we can not exclude other configurations of magnetic and orbital order, the 
ordered state certainly has large unit cell structure with a nontrivial periodicity concerning 
spin and orbital pattern. It would be desired to examine the spin-orbital order proposed in 
this paper in careful experimental studies. 

Since the Mott gap is very small, it would be interesting to apply pressure, because the 
gap may vanish and be metallized under reasonable and accessible range. Although the Mott 
transition through the bandwidth control route can be observed in many organic conductors, 
there exist not many examples in inorganic systems including transition metal compounds. 
If a large single crystal may be grown, it will greatly serve in clarifying the nature of the 
Mott transition since sample size of the organic conductors is in most cases too small to allow 
various experimental studies such as the neutron scattering. A large sample may serve in 
clarifying the nature of Mott transition in this route, where unusual criticality is proposed. 38 
Many of the organic conductors have superconducting phases near the Mott transition while 
the superconductivity is absent in metals close to the bandwith control Mott transition for the 
transition metal compounds such as V2O3 under pressure and LaTiC>3 with oxygen deficiency. 
Studies of metallic phases expected for Sr2V04 under pressure will contribute in clarifying 
the mechanism of superconductivity observed near the Mott transition including the cuprate 
superconductors. Another intriguing point is large orbital degeneracy and fluctuations possibly 
expected in the metallic phase of this material, which may have different aspect from other 
cases. In fact, our result for the optical conductivity and the comparison with the experimental 
result show that, in the insulating side, the insulating gap is selectively reduced for electrons on 
the xy orbital while the gap of electrons on the zx and yz orbitals remain. This is reminiscent 
of the orbital selective Mott transition recently proposed for (Sr,Ca)2RuC>4. 
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